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The turbulent dynamo effect, which describes the generation of magnetic fields in astrophysical 
objects, is described by the dynamo equation. This, in the kinematic (linear) approximation gives 
an unbounded exponential growth of the long wavelength part of the magnetic fields. Here we, in 
a systematic diagrammatic, perturbation theory, show how non-linear effects suppress the linear 
instability and bring down the growth rate to zero in the large time limit. We work with diffferent 
background velicity spectrum and initial magnetic field correlations. Our results indicate the robust- 
ness and very general nature of dynamo growth: It is qualitatively independent of the background 
velicty and intial magnetic field spectra. We also argue that our results can be justified within the 
bX). framework of the first order smoothing approximation, as applicable for the full non-linear problem. 

^ , We discuss our results from the view points of renormalisation group analysis. 
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I. INTRODUCTION 
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. ^ . Magnetic fields are ubiquitous. All astrophysical objects are known to have magnetic fields of different magni- 
tudes, e.g., 1 gauss at the stellar scale to 10~^ gauss at the galactic scale [Q. The origin of this field {primordial field) 
is not very clear - there are several competing theories which attempt to describe this However, a finite magnetic 
field in any physical system undergoes a temporal decay due to the finite conductivities of the medium. So, for steady 
magnetic fields to remain in astrophysical bodies, there has to be regeneration of the magnetic fields which takes place 
in the form of dynamo process . Typically astrophysical bodies are thought to have fast dynamo operating within 
themselves, resulting into exponential growth of the magnetic fields. This mechanism requires a turbulent velocity 
background [Q| [though non-turbulent velocity fields too can make a seed (initial) magnetic field to grow (for details 
see y]) we will not consider such cases here]. Since the dynamo equations, in the linear approximation (see below) 
gives unbounded exponentially growing solutions for the long wavelength (large scale) part of the magnetic fields, it is 
linearly unstable in the low wavenumber limit. An intriguing question, that arises very naturally is, what happens at 
a later time, i.e., whether magnetic fields continue to grow even after a long time. This obviously does not happen as 
we do not see ever growing magnetic fields in the core of the earth or in the sun. For example, geomagnetic fields (~ 
1 gauss) are known to be stable for about 10^ years 0]. Secondly, continuously growing magnetic fields violate energy 
conservation. In other words, if the dynamo equation correctly describes the problem, then its physically realisable 
^ solutions must not be linealy unstable in the long time limit. So there must be a counter mechanism to stabilise it, 
, which we investigate here. 

^ ' There have been numerous works in this field in the past by many groups. For examples Pouquet, Frisch and Leorat 
5^ 1 eddy damped quasi-normal Markovian approximation studied the connections between the dynamo process 
and the inverse cascade of magnetic and kinetic energies. Moffatt ||^, by linearising the equations of motion of three- 
dimensional (3(i) magnetohydrodynamics (MHD) examined the back reactions due to the Lorentz force for magnetic 
Prandtl number Pm ^ 1. Vainshtein and Cattaneo Q discussed several nonlinear restrictions on the generations of 
magnetic fields. Field et al discussed nonlinear a-effects within a two-scale approach. Rogachevskii and Kleeorin 
1^ studied the effects of an anisotropic background turbulence on the dynamo process. Brandenburg examined non- 
linear a-effects in numerical simulation of helical MHD turbulence [^. He particularly examined the dependences 
of dynamo growth and the saturation field on P^- Bhattacharjee and Yuan pO[ | studied the problem in a two-scale 
approach by linearising the equations of motion. However these issues are not yet fully closed. We examine the 
following questions: i)instead of a two-scale approach (which is rather adhoc) whether we can employ a diagrammatic 
perturbation theory, which has been highly successful in the context of critical dynamics [ pT| , driven systems etc., 
can be easily extended to higher orders in perturbation expansion and provides natural connections with standard 
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renormalisation group framework, and ii)if a turbulent background |^ is essential for the dynamo mechanism. To 
put it differently we ask if dynamo process can take place with velocity fields with arbitrary statistics. We explicitly 
demonstrate that the nonlinear feedback of the magnetic fields on the velocity fields in the form of the Lorentz force 
stabilises the instability. We show it for a very general velocity and initial magnetic field correlations - thus our results 
demonstrate the very generality of the dynamo process. The plan of the rest of the paper is as follows: In Section || 
we discuss the general dynamo mechanism within the standard linear approximation. In Sec. Ill B| we show that one 
needs to go beyond the linear approximation, i.e., include the non- linear effects to see eventual saturation of magnetic 
field growth. In Sec. |^ wc conclude. 



II. DYNAMO GROWTH: THE LINEAR APPROXIMATION 

In the kinematic approximation i.e., in the early time when magnetic energy is much smaller than the kinetic 

energy (/ u^d^r >> J b^d^r), where u(r,t) and b(r,i) are the velocity and magnetic fields respectively, the Lorentz 
force term of the Navier Stokes equation is neglected. In that weak magnetic field limit, which is reasonable at an 
early time, the time evolution problem of the magnetic fields is a linear problem as the Induction equation is 
linear in magnetic fields b: 

^ = Vx(uxb) + AiV2b, (1) 

where fj, is the magnetic viscosity. The velocity field is governed by the Navier-Stokes equation (dropping the 
Lorentz force) 

^ + u.Vu = -— + vV^u + {. (2) 
ot p 

Here v is the fluid viscosity, f an external forcing function, p the pressure and p the density of the fluid. We take f 
to be a zero mean, Gaussian stochastic force with a specified variance (see below). 

In a two-scale ^ approach one can then write an effective equation for B, the long- wavelength part of the magnetic 
fields §]: 

— = V X (U X B) + V X E + ^V^B, (3) 
ot 

where the Electromotive force E = (v x b). U is the large scale component of the velocity field u. An Operator 
Product Expansion (OPE) is shown to hold jl^ which provides a gradient expansion in terms of B for the product 
E = (u X b) [| 

dB ■ 

E, = a,j Bj + (3ijk-g^ + ■■■■ (4) 

For homogenous and isotropic flows (a^ — aSij) Eq.(^) gives, 

— = V X (U X B) + aV X B + /iV^B, (5) 
ot 

which is the standard turbulent dynamo equation. Here /i now is the effective magnetic viscosity which includes 
turbulent diffusion, represented by /^yj, in Eq.(^. a depends upon the statistics of the velocity field (or, equivalently, 
the correlations of f). Retaining only the a -term and dropping all others from the RHS of Eq.(||), the equations for 
the cartesian components of B become (we neglect the dissipative terms proportional to as we are interested only 
in the long wavelength properties) 



^By a turbulent background we do not mean any kind of fluctuating state but a fluctuating state with K41 spectrum for the 
kinetic and magnetic energies and cascades of appropriate quantities; if there is no mean magnetic field then the energy spectra 
is expected to be K41-type - see Ref. Q . 
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B,(k,t)\ / ~k, ky \ /B,{k,t)\ 

By(k,t) \ =ia\ k, -fc, By{k,t) 
B,(k,t)J \-ky k^ J \B,{k,t)J 



The eigenvalues of the matrix is A = ±zfc, 0. Thus depending on the sign of the product afc, one mode grows and 
the other decays. The third mode stays constant in time. Since growth rate is proportional to |fc| and dissipation 
proportional to fc^, large scale fields continue to grow leading to long wavelength instability. Thus in the long time 
limit effectively only the growing mode remains. Since the cartesian components of B are just linear combinations 
of the eigenmodes of the matrix above, they also grow exponentially. Growth rate a is a pseudo-scalar quantity, i.e., 
under parity transformation r —f — r, a —>■ —a [p]p^]. Since a depends upond the statistical properties of the velocity 
field, its statistics should not be parity invariant. This can happen in a rotating frame, where the angular velocity 
explicitly breaks parity. 



III. FORMULATION OF THE DYNAMO PROBLEM IN A ROTATING FRAME 

The NS and the Induction equation in the inertial (lab) frame in (k, t) space become 

^^^^ + ^Pvpi^) E "^ (q- - q, t) = ^Pyp(k) E t)bpi^ - q, i) + ^V^i.. + /.(k, t), (6) 



db,{k,t) 
dt 



Here, Ui{k,t) and 6j(k, t) are the fourier transforms of Ui{r,t) and bi(r,t) respectively, Pijp{k) = Pij{k)kp + 
Pip(k)kj, Pijp{k) = Pij{k)kp — Pip{k)kj, Pij is the projection operator, which appears due to the divergence-free 
conditions on the velocity and magnetic fields. The Eqs.(^ and are to be supplimented by appropriate correla- 
tions for fi and intial conditions on bi. We choose /i(k, t) and 6i(k, t = 0) to be zero mean Gaussian distributed with 
the following variances: 

{Mk,t)f,i-k,0)} ^ P,,D,ik)6{t), (8) 

{b,{k,t)bj{-k,0)) = P,,D2{k)6{t). (9) 

Di and D2 are some functions of k (to be specified later). 

In a rotating frame with a rotation velocity $1 = ilz the Eqs.(^ and (|^) take the form 

g^' + 2(fl X u), + ]^P^,p[k) E u,{ci, t)up{k - q, t) = ]^P^,p{k) ^ (q, t)bp{k - q, i) + yV^u^ + /,(k, t), (10) 



^^'^^'^^ ■ (nxb), = P„p(k)VMj(q,<)fcp(k-q,<)+/iV2&,. (11) 
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r2 X u is the coriolis force. The centrifugal force x (Jl x r) is put in as a part of the effective pressure=p - 
which does not contribute to the dynamics for incompressible flows. The correlations Eqs.(^) and (||) do not change 
in the rotating frame. The bare propagators Gu and Gb of Ui and bi are 
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such that u = Gu f and b(<) = Gb h{t — 0) where 
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We get for the correlation function (Mi(k, k, — w)) = (uu-'") ( u-^ is the transpose of u), when i ^ j there are 

terms proportinal to 17 (for i — j, there are no 0{Q) terms). Under parity transformation these i ^ j terms change 
sign. Thus the i ^ j terms have both even parity and odd parity parts. Similarly in (bibj), fori ^ j have both odd 
and even parity parts. A direct consequence of these parity breaking parts is that fluid helicity (= J^(u.(V x u))) is 
non-zero: /^(u.(V x u)) oc fl. In the same way /^(b.(V x b)) is non-zero and proportional to fl. Both J^(u.(V x u)) 
and /^(b.(V x b)) have same signs - a fact of great importance for the results discussed here. Notice that G^'^ are 
different from G^'^ yy - this is just the consequence of the fact that Q. distinguishes the z-direction from others, making 
the system anisotropic. However for frequencies w >> il or length scales >> VL isotropy is restored. In that regime 
the role of the global rotation is only to introduce parity breaking contributions proportional to to (uiUj) and {bibj) 
for i j. These can be modeled by introducing parity breaking parts in Eqs.(^) and (|^) 

(/,(k,i)/,(-k,0)) = P,,Di{k)d{t) + 2ie,,pkpDiik)S{t), (12) 
(&,(k,t)6j(-k,0)) - P,jD2{k)5{t) + 2ie,,pkpD2ik)5{t), (13) 

in conjunction with the inertial frame Eqs.(^ and (Q). The parity breaking parts in the noise correlations or intial 
conditions ensure that /^(u.(V x u)) and J^(b.(V x b)) are non-zero as is the case with Eqs.([lO|) and Eqs.([rT|) along 
with Eqs.(^) and (^). What is the relative sign between Di and -D2? Since J'^(u.(V x u)) and J^(b.(V x b)) are 
proportional to Di and D2 respectively, and they have same signs, Di and D2 must have same signs. 



A. Calculation of a in the kinematic approximation 



In the kinematic approximation neglecting the Lorentz force term of the Navier-Stokes equation, the time evolution 
of the magnetic fields follows from the linear Induction Equation (|l|). We assume, for convenience of calculaions, 
that the velocity field (u) statistics has reached a steady state. This is acceptable as long as the loss due to the 
transfer of kinetic energy by the dynamo process is compensated by the external drive. In the kinematic (i.e., linear) 
approximation, we work with the Eqs.(^ (without the Lorentz force) and @j. We choose /z(k, t) to be a zero-mean, 
Gaussian random field with correlations 

(/i(k,i)/„,(k,0)) = 2PirnDiik)S{t) + 2ieimnDi{k)kr,Sit). (14) 

Our intial conditions for the magnetic fields are 

{ba{K t = 0)6^(-k, i = 0)) - Pa02D2{k) + 2ie^p^k^b2{k), (15) 

Since we are interested to investigate the dynamo process with arbitrary statistics for the velocity and magnetic fields 
we work with arbitrary Di{k), Di{k), D2{k) and D2{k)- For K41 spectra, we require |l^ Di{k) — Dik^^, Di{k) = 
Dik~'^ , D2{k) = D2k^^/^ and D2{k) = k^^/^. These choices ensure that under spatial rescaling x /x, v,b ^ 
b} which is the Kolmogorov scaling. Starting with such a correlations ensures that only the amplitudes of the 
magnetic field correlations grow in time due to the dynamo effects but the inertial range scale dependence does not 
change. However, this may not be the case always. In general, not only the amplitude grows in time due to dynamo 
actions, the scale dependence too can be anything at i = and evolve in time. Note that both the force correlations 
in the Eq.(^ and the intial conditions on Eq.(|^) have parts that are parity breaking, in conformity with our previous 
discussions. We now calculate the a-term (at the tree level) in the kinematic approximation (which we call the 'direct' 
term - responsible for growth) below (see Fig. la): 

((u X b)^)D = ( / e^/3T,w^(q,t)6^(k - q,<)} ( / eQ^^w^(q, t)e^5Ai(k - q)5W^(qi, ti)6T(k - q - qi, ii)G'o(k - q, i - ti)) 

(16) 

which gives the a-term: 
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a_DBa(k,t) 



1 exp(- 2*1/92) 

q^iv + li) q^{iy-fi) 



giving an 



1 f 2_2(£l_ 

3 l'(v+p) Jq (u+fj,)q-' 



2S3 



for large t. Thus self-consistently, a_D = — ^ 2 ^j | ^ 



_DiM_ 



8|-(k+p)'PT' 



(17) 



The suffix D 



refers to growth or the direct term, as opposed to feedback which we discuss in the next Sec. [II B . The growth term is 
proportional to and diffusive decay proportional to k^. So large scale components grow and small scale components 
decay. 



O 



Di 



O 



/ - q 



k-q k k-q k 

FIG. 1. Tree level diagrams for < u(q) x b(k — q) >. (a) A solid line indicates a bare magnetic field response function, 
a broken line indicates a bare velocity response funtion, a 'o' joined by two broken lines indicates a bare velocity correlation 
function (proportional to -Di), a wavy line indicates a magnetic field, a solid triangle indicates a ub vertex. This contributes of 
an (b) A solid line indicates a bare magnetic field response function, a broken line indicates a bare velocity response funtion, 
a 'o' joined by two broken lines indicates a bare magnetic field correlation function (proportional to D2), a wavy line indicates 
a magnetic field, a solid triangle indicates a ub vertex. This contributes to ap- 



B. Suppression of growth rate: Nonlinear feedback 

When the magnetic fields become strong neglecting the feedback of the magnetic fields in the form of the Lorentz 
force is no longer justified. So we work with the full Eqs.(||) and (0). We here follow a diagrammatic perturbation 
approach. In presence of the Lorentz force there is an additional contribution to a (Fig. lb). 



((u X bj)i^} = ( / e.,jpUj{q,t)bp{k - q,t)) (18) 

Jq 

P,ran{q)Gl{q, t - ti) 5„ ((ji , i i )6„ (q - qiM)Gl{k - q, t)hp{k - q,t = 0)) (19) 
which gives {F refers to feedback) 

aB,ik,t)^^e,,p f P,-„.„(g)c2"°H*-Wt5„(k^ ,) -2^D,{q)e q 
Jq 2aD\q\~2m 



(20) 



This on simplification gives 



3 15 L a_D|q| - 2/^(j-^ 



where D2{q,t) — exp[2aD\q\t — 2iiq^t]D2{q) is a growing function of time for small wavenumbers. Thus ap grows in 
time. 

Thus, at a late time t, when the non-linear feedback on the velocity field due to the Lorentz force is nolonger 
negligible, i.e., for a finite ajj and ap we find self-consistently, 

^^^[^ EM (21) 

^ 3 J (2^)3t.[|(a^-t-a^)g|-(^. + ;,)q2]' ^ } 

_ 2^3 4 f d\ D2{q,t) 

3 15 7 {27rr\iaD + aF)q\~2tiq^- ^ ' 
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-^ = ^^J 7^ .u.._ , AV„, o..„J - (25) 



Thus net growth rate oc \{aD + aF)k\ for the mode i?i(k, t). 

Let us now consider various k dependences of Di{k) and D2{k). For the case when the background velocity field is 
driven by the Navier-Stokes equation with a conserved noise (thermal noise) one requires to have Di{k) = Dik'^, Di ~ 
Dik, giving (wi(k, k, t)} = constant. If we assume similar fc-dependences for (6i(k, 0)6i(— k, 0)) then we require 

1^2 (fc) ^ constant and D2{k) = These choices give 

^^^[^ Diq 1 .24^ 

3 J (2^)3i.[|(a^+a;^)g|-(^. + ^)q2J' ^ > 

2S3± f ^ D2it) 

3 15 J {2tt)^ q[\{aD + aF)q\ ~ 2fiq^ 

which remain finite even as the system size diverges. 

Fully developed turbulent state, characterised by K41 energy spectra, is generated by Di{k) ^ k^^ and -Di(fc) — 
Dik~'^. In addtion if we assume that the initial magnetic field correlations also have fc41 scaling then D2{k) ^ k^^/'^ 
and D2(k) = D2k~^/^ . If one starts with a K41-type initial correlations for the magnetic fields, then at a later time 
the scale dependence for the magnetic field correlations are likely to remain same; only the amplitudes grow. Notice 
that the spectra diverge as fc ^ 0, i.e., as the system size diverges. This is a typical characteristic of fully developed 
turbulence. For such a system self-consistently, 

^^^Ssf^ Dig-* 

3 J i27rr ,y[\iaD + aF)q\ ~ i:^ + fi)q'y ^ ' 

_ 2^3 4 f d^q D2it)q~'^' 
^ 3 15 7 {27rr [\{aD + aF)q\ - 2nq^Y ^ 



The notable difference between the expressions Eqs.(g5|) and Eqs.(|27|) for the a coefficients is that the a coefficients 
diverge with the system size when the energy spectra are singular in the infra red limit (for fully developed turbulence). 

In general, at early times (small uf), ap increases exponentially in time. For t ^ ^Ina^), ao and ap are 
comparable. The growth rate of q;f comes down and a^i decays. Since au and have different signs, |(a£)+aF)| — > 
as t — > large. Thus the net growth rate comes down to zero. Hence, Eqs.(p5|) and Eqs.(^) suggest that the early time 
growth and late time saturation of magnetic fields take place for different kind of background velicity correlations 
and initial magnetic field correlations. Thus dynamo instability and its saturation are rather intrinsic properties of 
the 3(iMHD equations in a rotating frame. Our results also suggest that these processes may take place for varying 
magnetic Prandtl number ^/v. The above analysis crucially depends on the fact that ao and ao have opposite signs, 
which, in turn, imply that Di and D2 have same signs. We have already seen that in a physically realisable situation 
where parity is broken entirely due to the global rotation, Di and D2 indeed have the same sign. Our results also 
suggest that the net growth rate is proportional to the difference between the fluid and the magnetic helicities at any 
time t. 

How can we understand our results in a simple way? To do that we resort to first order smoothing approximation 
In the kinematic limit, in this smoothing approximation to calculate (v x b) one considers only the Induction 
equation as v is supposed to be given. However when one goes beyond the kinematic approximation, one has to consider 
the Navier-Stokes equation as well. Thus in the first-order smoothing approximation one writes the equations for the 
fluctuations v and b as (to the first order) 

^ « V X (v X B) + V X (V X b), (28) 
ot 

and 

|^«... + (B.V)b, (29) 

where the ellipsis refer to all other terms in the Navier-Stokes equation and B and V are the large scale {mean field) 
part of the velocity and magnetic fields . With these we can write 

d d 

(v X b}j = {eijpVjbp) {<^ijpVjB.m- — Vp) + {eijpbpBmj. — hj) = UimBm + ■■ ■ (30) 

OXm CfXrn 

Here the ellipsis refer to non-a terms in the expansion of (v x b) (see Eq.(4)). Thus for isotropic situations a = 
(v.(V x v)) + (b.(V X b))]. Thus a is proportional to the difference in the fluid and magnetic helicities, a result we 
have already obtained through a more detailed calculation above. In our model fluid helicity is statistically constant 
in time, but magnetic helicity grows in time and hence a — > in the long time limit. 
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IV. SUMMARY 



Thus in conclusions we have shown how the initial exponential growth of the magnetic fields in a turbulent dynamo 
can be arrested by the action of the magnetic fields on the velocity fields. Our mechanism required that the parity 
breaking parts of the velocity and magnetic field variances must have the same sign, which must be the case in any 
physical system. It is also worth noting the role of the symmetries of the velocity and magnetic field correlation 
functions. The antisymetric part, which is also an irreducible part of the (viVj) tensor is responsible for the growth. 
Again the antisymetric (irreducible) part, part of the (bibj) tensor is responsible for stabilisation. Even though our 
explicit calculations were done by using simple initial conditions for the calculational convenience, the results that we 
obtain are general enough and it is obvious that the feedback mechanism is independent of the details of the intial 
conditions. Thus our results should be valid for more realistic initial conditions also. We have also demonstrated 
that our results can be easily understood within the first-order smoothing approximation. From the point of view 
of nonlinear systems, our results can be interpreted as 'non-linear stabilisation of linear instabilities', qualitatively 
similar to the well-known example of the Kuramato-Shivanisky (KS) equation in one and two dimensions. This is 
linearly unstable for small wavenumbers, but the nonlinear term stabilises it. In fact, the long wavelength properties 
of the KS equation is same as that of the stochastically driven Kardar-Parisi-Zhang equation (KPZ) ||l9| in one 
and two dimensions. However we must add words of caution while making the comparison between the KS-KPZ 
problem and the dynamo problem. In the KS-KPZ case the KS equation is not stochastically driven. The long 
wavelength instability serves as drive. However, in the present case of dynamo, the velocity field (or the NS equation) 
is stochastically driven. Recently it has been shown that the stochastically and determinitically driven NS equations 
belong to the same multiscaling universality class ||2l[] (in the inertial frames). Even though the same has not been 
shown for 3dMHD, it is probably true there too. So with some confidence we can draw the analogy between the 
KS-KPZ problem and the dynamo problem which we discussed here. An important issue is still however left open. 
In fully developed SdMHD, in the steady state, correlation and response functions exhibit dynamical scaling and the 
dynamic exponent z — 2/3 |2^ , ^ , which means renormalise dissipations (kinetic as well as magnetic) diverge ~ k~'^/^ 
for a wavenumber k belonging to the inertial range. Even for decaying MHD with initial K41-type correlations this 
turns out to be true where equal time correlations exhibit dynamical scaling with z = 2/3. The question is, 
what it is in the intial transient of dynamo growth? K41-type of spectra suggest the roughness exponent for both the 
velocity and the magentic is 1/3 This, together with Galilean invariance for the 3(iMHD equations give 

z = 2/3. So, unless a coefficients too pick up divergent corrections dynamo growth will be subdominant to dissipative 
decay in the long wavelength limit. However, as our results [Eqs.(|^] suggest, the alpha coefficients diverge in the 
long wavelength limit, indicating that they pick up fc-dependent singular corrections. Simple minded calculations 
[0 suggest a ~ fc^^/"^. If this is really the case then neither growth nor dissipative decay dominate in the inertial 
range; the sign of (a — /i) determines it. However this kind of RG calculations suffer from few technical problems 
[ pTf . Thus to settle it conclusively one requires more sophesticated technique and/or numerical simulations. This 
issue of divergent effective viscosities in the inertial range assumes importance as it may help to overcome some of the 
non-linear restrictions as discussed by Vainshtein and Cattaneo . A system of magnetohydrodynamic turbulence in 
a rotating frame, after the saturation time (i.e., after which there is no net dynamo growth) belongs to the universality 
class of usual three-dimensional magnetohydrodynamic turbulence in a laboratory. This can be seen easily as both in 
the lab and rotating frames, the roughness and the dynamic exponents can be calculated exactly by using the Galilean 
invariance and noise-nonrenormalisation conditions p^j2^ . However, the same cannot be immediately said about the 
multiscaling exponents of the higher order structure functions. Further investigation is required in this direction. It 
will also be very interesting to find out the detailed quantitative dependences of a on the magnetic Prandtl number 
(Pm) in view of the recent results that Pm is connected to other dimensionless numbers like the ratio of the total cross 
helicity to the kinetic energy in the steady state [p5| . 
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